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I. INTRODUCTION 

The purpose of this paper is (i) to revisit the formulation of the 
relativistic Sunyaev-Zel'dovich (rSZ) distortion of the cosmic 
microwave background (CMB) due to a population of non- 
thermal high energy electrons in clusters of galaxies, and (ii) 
to provide a unified and more comprehensive framework to 
compute the rSZ effect. 

Inverse Compton scattering of non-relativistic thermal elec- 
trons on cosmic microwave background (CMB) photons has 
been widely investigated in the literature since the 1960- 
1970's, notably in the context of the seminal works by Sun- 
yaev and Zel'dovich on the SZ effect 11112], and in early stud- 
ies on relativistic electron cosmic rays in high energy astro- 
physics (e.g. iH 0, Htl). The SZ effect generated by a pop- 
ulation of hot (still non-relativistic) electrons has also been 
addressed in detail since the 1980's (see e.g. 0. IE & 
[M [H [H [il B), while the SZ effect originating from 
relativistic electron cosmic rays has been discussed in e.g. 
Refs. H H [H [Htl. So far, these works on hot thermal 
or relativistic non-thermal electrons have used different ap- 
proaches. 

For example, in Refs. Hol [TH [ll [Tl [Till , the authors have 
computed the relativistic corrections to the standard thermal 
SZ predictions due to the hot trail of intracluster thermal elec- 
trons (including the effect from multiple scattering). They 
used a Fokker-Planck expansion of the collisional Boltzmann 
equation in terms of the electron temperature-to-mass ratio 
r/m < 0.1, in which the full squared matrix amplitude of 
Compton scattering (that is given in any textbook, e.g. (O) 
was implemented for the relevant kinematics. Such an ap- 
proach therefore not only encompasses but also extends the 
usual calculation of the standard thermal SZ effect, thus mak- 
ing this formalism very attractive. Nonetheless, it has not 
been applied yet to the case where there are two populations 
of electrons (one dominant semi-relativistic population and a 
subdominant relativistic population). 
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The authors of e.g. Refs. ||lH[i[il[i3,[il have considered 
the SZ contributions coming from a population of thermal 
electrons plus an extra source of non-thermal electron cosmic 
rays, as we want to do. However, they used a more empir- 
ical approach based on a classical radiative transfer method 
proposed by Chandrasekhar ll20ll . Indeed, they computed the 
spectral shift in the CMB photon intensity due to Compton 
scattering processes with electrons, scattering off some pho- 
ton frequencies (from energy E to £') and scattering in some 
others (from energy E' to E). The squared matrix amplitude 
(that is the probability of frequency change) was expressed 
in the electron rest frame, leading to an expression close to 
its associated non-relativistic limit, such as the one originally 
used by Chandrasekhar. Multiple scatterings were eventually 
implemented by further using Chandrasekhar's radiative ap- 
proach. 

In Sect. |II] we demonstrate that Chandrasekhar's radiative 
transfer method and the Boltzmann-like approach mentioned 
above are in fact formally equivalent in the single scatter- 
ing approximation. They thus give the same results, which 
we show by using the full relativistic squared matrix ampli- 
tude instead of the photon frequency redistribution probabil- 
ity function derived by Wright [63. This therefore sketches 
a unified formalism for the computation of both hot thermal 
and relativistic non-thermal SZ effects, as subsequently found 
in ll2lll . For illustration, we derive the relativistic photon en- 
ergy redististribution function in the CMB reference frame 
(therefore different from Wright's P{s,^) function) that can 
be used to compute the relativistic SZ effect generated by an 
additional population of high energy electron cosmic rays in 
a cluster. For the sake of completenes, we revisit in Sect. [Ill] 
the full relativistic phase space and squared amplitude of the 
inverse Compton process in the relativistic limit. Our conclu- 
sions are drawn in Sect. |IV] Numerical applications of this 
framework will be given in a companion paper ll22ll . 



II. EQUIVALENCE BETWEEN RADIATIVE TRANSFER 
AND THE BOLTZMANN EQUATION 



The purpose of this section is to demonstrate the equivalence 
between the approach based on the Boltzmann equation and 
Chandrasekhar's radiative transfer method to compute the rSZ 
effect, in the single scattering approximation. 
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A. Scattered intensities in the radiative transfer approach 



B. Boltzmann-Hke formalism 



Chandrasekhar's radiative transfer method has been widely 
used in the hterature, notably to compute the SZ effect due 
to non-thermal electrons in clusters of galaxies (for a review, 
see lfl6ll ). The principle of this approach is recalled in the ap- 
pendix (see Sect. |A]i, as well as the precise derivation of the 
CMB intensity deviation in the relativistic regime and in the 
single scattering approximation (see Sect. |B]i. Here we sum- 
marize the main results coming out of this approach. 

The spectral shift in the CMB intensity due to scatterings with 
relativistic electrons can be expressed as the difference be- 
tween two contributions. The first one is associated with the 
photons scattered out (/^^""'(E^)) from the CMB spectrum, 
from energy E/^ to energy ^ Ei^, and the second is associ- 
ated with the photons scattered in (/^^'"(iii)) from energy E'f, 
to energy Ek (after averaging over the direction of the incom- 
ing photons): 



AlyiEk) 



(1) 



The expressions of the in and out intensities can then be cal- 
culated. They are given in Sect. |B] (see Eqs. IB2l and IB3l l as 
integrals over the line of sight dl: 
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where ly is the intensity associated with the photon black- 
body spectrum (before interaction with the electrons) and fe 
is the energy distribution function characterizing the relativis- 
tic electrons (of energy Ep). The term |9l/ p is the squared 
matrix amplitude associated with the process {p,k — * p',^') 

for the scattered-out part, and |sv/p describes the process 
{p,k' — > p',k) for the scattered in part. The terms t and t 
feature the photon energy transfer for both processes, respec- 
tively, and are defined in Eqs. (l25IB3l l. /j characterizes the 
angle between the incoming electron of momentum p and the 
photon of momentum k (see Eq.|27t. 

As we will show in the next section, the photon frequency 
redistribution probability function considered by Wright in 
1979 lit], which has been continuously used afterwards for 
rSZ predictions H [H [H [H, will differ from our analy- 
sis in that it was derived in the electron rest frame, but both 
are consistent ll2lll . Anyway, the equivalence between the 
radiative transfer and covariant approaches — in the relativis- 
tic regime and single scattering approximation — can be for- 
mally demonstrated without the need for replacing the squared 
matrix amplitude by its expression. 



1. Covariant collisional Boltzmann equation 

The theoretical framework developed in ifTol [TTl [l3l [Till to 
treat the semi-relativistic corrections of the thermal SZ is the 
covariant collisional Boltzmann equation (see e.g. ll23ll24ll25i 
|26ll ), which expresses the momentum transfer between differ- 
ent interacting species. 

When focusing on photons scattered by a gas of electrons in 
a covariant scheme, the spatial convective currents are irrele- 
vant, and the time evolution of the photon phase space is given 
by the covariant collisional Boltzmann equation as: 
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The collision integral is defined by: 
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where fy and are the photon and electron distributions, re- 
spectively, and where we have neglected the Pauli blocking 
and stimulated emission factors (even if the latter were in- 
cluded, they would automatically disappear as soon as treat- 
ing relativistic electrons, since in this case fe{Ep) ~ fe{Epi), 
as we will show in the following). This equation features the 
Lorentz invariant phase space c/LIPS and the squared matrix 
amplitude |fAf P associated with the process [p^k ^ p',k'). 
These will be made explicit in Sect.|IIIl Note that this formal- 
ism rests on the microreversibility of the process. 
Let us now discuss the meaning of the electron distribution 
featured above. As we have already stressed, we are consid- 
ering two different populations of electrons, one thermal and 
the other non-thermal. The former is assumed dominant, but 
we are interested in the subdominant population only. For- 
mally, the above electron distribution characterizes the sum 
of these two populations, that is fe = /"^* +/*■ Hence, the 
question of whether we can separate their respective contri- 
butions arises. When integrating over the electron momenta, 
we see that at low energy, say below 1 MeV typically, only 
the thermal part /* will be relevant in the calculation, while 
at higher energy, only the non-thermal part will contribute. 
Since we restrict ourselves to the single scattering limit, we 
can therefore safely separate those two contributions since 
none of them will affect the other one. Consequently, in the 
following, we will only deal with the relativistic non-thermal 
part, and will set fe = /"^* for convenience. Note that in the 
multiple scattering case, this would not hold anymore, since 
there would be a mixing of the SZ effects coming from both 
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populations of electrons: for instance, a photon could first be 
scattered by a thermal electron, and then by a non-thermal 
one. On that account, we would have to treat both popula- 
tions in a self-consistent manner. Nevertheless, such a com- 
plex analysis is beyond the scope of this paper, and probably 
too sophisticated regarding the small contribution to the SZ 
effect that is expected from relativistic electrons. 
The collisional Boltzmann equation written above describes 
the time evolution of the photon phase space distribution due 
to elastic collisions with electrons. The electron distribu- 
tion has therefore to be known at any time. In the following 
(Sect. Ill B 21 ). we will demonstrate that if one assumes that 
the electron phase space density does not change significantly 
over one (or even several) characteristic travel time of a pho- 
ton through a cluster, then one can use this formalism to treat 
the SZ effect originating from non-thermal electrons in the 
single scattering approximation. The above hypothesis is jus- 
tified since, if we take a typical scale of 10 Mpc for a cluster, 
then the mean time for a photon to cross it is ~ 30 Myr while 
the typical timescale for the relativistic electron energy loss by 
inverse Compton scattering on CMB photons is ~ 300 Myr 



2. Covariant formalism applied to the rSZ ejfect 

Let us consider a cluster containing a tiny population of elec- 
trons, and track the photon distribution. The integration over 
the time, in Eq. (|2]), translates to an integration over a geodesic 
(dl = c dt): 
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If we consider the electron density to be small enough, we 
can assume that a CMB photon has a very small probability 
to experience Compton scattering with electrons along its path 
through a cluster. In this case, the single scattering approxima- 
tion is fully relevant and gives; / dfy = f^^^'^' — fy^^°^^, where 
ybefore jj^g original black-body spectrum distribution before 
the CMB photons enter into the cluster, namely fy, and /y*^'" 
is the photon distribution after travelling through the cluster 
and experiencing one scattering. An observer looking in the 
direction of the cluster will therefore observe fy^^^'^ ~ fy^ fy- 
Expressing now this equation in terms of the photon intensity. 



t\\?A\fiIy{E)=E^fy{E), we obtain: 
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where is the intensity of the pure black-body spectrum of 
the CMB before entering in the cluster. The first term in the 
sum under the integral characterizes the frequencies shifted 
out from Ek to Ef^i, while the second term characterizes the 



frequencies shifted in from E^-i to Ek, both due to Compton 
scattering. Note that the electron distribution /g is not ubiq- 
uitous anymore, and must now also carry information of the 
spatial distribution, which will be relevant for the line of sight 
integral along dl. 

As we stated previously, if the number density of the rela- 
tivistic electron population is very small, then the interactions 
between the CMB photons and this population is well approx- 
imated by the single scattering approximation. This hypothe- 
sis further enables us to treat several populations of electrons 
(e.g. a dominant, thermal and/or semi-relativistic, population 
and a subdominant relativistic population) separately. 
Let us now specifically focus on the relativistic population. 
As we are considering electrons with typical energies Ep>l 
MeV and CMB photons with energies Ef^ < 10^^ eV, we can 
safely make use of the limit a = E^/Ep ^ in the above 
equation, which is fully valid in the whole range of the CMB 
photon energy distribution in the single scattering approxima- 
tion. Within this limit, due to energy-momentum conserva- 
tion in the Compton scattering, we will show in Sect. Hill that 
Epi Ep, such that we can replace fe{Epi) by fe{Ep) (see 
Eqs.[33]and[35]). This is the key point to understand why the 
Boltzmann framework can be linked to the radiative transfer 
approach for the SZ effect associated with non-thermal elec- 
trons. Most importantly, in the single scattering approxima- 
tion, the photon intensity ly under the integral of the right hand 
side can be replaced by that corresponding to the pure black- 



body spectrum /5. Defining My 



My{Ek) 




- ly, we obtain; 
fe{Ep,x) c/LIPS 



(6) 



where t = Ej^i /Ek, and where we have explicitly added a fac- 
tor of two to account for the sum over the two photon polar- 
izations (so that it should not appear in the definition of ly 
in the right hand side). This expression fully characterizes 
the intensity after removing the contribution of photons scat- 
tered out (from Ek to Ej^i) and after adding that of photons 
scattered in (from to Ej^), exactly like a la Wright. The 
term proportional to Iy{Ek) is the scattered out contribution, 
and the term proportional to Iy{Ei-i) is the scattered in contri- 
bution. Hence Chandrasekhar's formalism, when applied to 
the relativistic SZ effect, as recalled in the previous section, is 
strictly equivalent to using the Boltzmann equation (in the sin- 
gle scattering approximation) — detailed expressions related 
to Chandrasekhar's formalism are derived in the appendix, see 
Sect.|B] We can rewrite the above equation in terms of the dif- 
ferential Compton cross section, after performing the integrals 
over d^p' and dEi^r. 
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where Piei is the relative velocity between the incoming parti- 
cles (see Eq.lSTTi. where t/ii^/ is the solid angle associated with 
photons of energy Ej^i ^ Ek, and where the differential Comp- 
ton cross section da is defined in Eq. ( |52] |. We will sketch the 
full angular dependence in the next subsection. 



C. Equivalence between the radiative transfer and covariant 
approaches 

Eqs. ^ or Q contain two contributions, which can conve- 
niently be expressed as: 



AIy{Ek)=l';{Ek)~I"''\Ek), 



(8) 



where the upper script out characterizes the part of the spec- 
trum which is shifted out from energy E/^ to other energies, 
while in stands for other frequencies which are shifted in to 
energy E/^ because of scattering. According to Eq. (|7]), the 
term Iy"^{Ek) is given by: 



di 



dO.^, p 



d^p 
da 
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where we recall that the factor of two in front of the right hand 
side accounts for the sum over the two photon polarizations 
(see Eq.|6]). This can be rewritten as: 



2xI^iEk), 



where we define: 
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We have used the definition of the differential cross section 
from Eq. (l53T l. and that of the relative velocity prei = (1 ^ 
P/j). Now we need to average over the solid angle dSl^ = 
d<^d/j associated with the incoming photons, so that the actual 
measured shifted intensity will be: 



\E,) 



2 T 

d£h^ 
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If we perform the integral over dQ.iti = dis/dp' appearing in the 
definition of T, we obtain the full relativistic (inverse) Comp- 
ton cross section given in Eq. (|54] |. which still depends on /j, 
Ek and Ep. Then, we can perform the average over the in- 
going photon solid angle dUf^ = 27tc//j. The result is found to 
depend on powers of ay^ = yEi^/m. Though for CMB photons 
ay^ is negligible up to ultra-relativistic electron energies, that 
is as long as y < 10^, we still give the result at next-to-leading 
order (see Eq.lSSTi: 



z = dl 
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where aj is the Thomson scattering cross section. Note that it 
is quite consistent to find a result in which the cross section is 
close to the Thomson cross section for the diffused out term, 
since we are in the limit yEj <C m. Indeed, this corresponds 
exactly to seeing a low energy photon in the electron energy 
rest frame: in this sense, as very well known, we are in the 
Thomson regime, while not dealing with a Thomson diffusion 
at all. If we had found other dependences on angles under 
the integral, it would have been more difficult to integrate the 
cross section out in the CMB rest frame, as it will be shown to 
be the case for the scattered in a bit further Note finally that a 
mere transformation of the above expression by re-expressing 
the zenith angles of the photons /j and /j' in the electron rest 
frame would have resulted in a more convenient form, as it is 
well explained in ll2lll . 

We recall that the spatial dependence of the electron popu- 
lation has to be included in the distribution fe- z can be, 
as displayed above, interpreted as the optical depth associ- 
ated with the electron population. To provide an expression 
closer to the standard lores, we can assume that the energy and 
the spatial distributions of electrons can be factorized, such 
that f(,{Ep,x) ~ ne{x)fe{Ep), where fe is the normalized mo- 
mentum distribution (/ d^p/ {2Tif'fe{Ep) = 1). This naturally 
gives / d^p/{2'!if'fe{Ep^x)~ne{x),io which one can stick in- 
stead of the previous trick, depending on one's convenience. 
We can therefore rewrite the previous equation as: 



(10) where: 
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(15) 



We thus recover the standard definition of the optical depth 
Tni for non-relativistic thermal electrons when neglecting the 
terms proportional to a and ay^. The factor 3C < 1 has to be 
considered as a relativistic correction to apply to Tm-. Finally, 
we remind that this optical depth will have to be averaged over 
the spatial resolution of the telescope for quantitative predic- 
tion purposes (see Eq.l23Tl. 

The scattered in term of Eq. (O, f^{Ek), is given, according 
to Eq. by: 



^ 2 dl 
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We recall that t = E^i jE^ is the ratio of the diffused to initial 
photon energies for the scattered-out process (and conversely 
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for the scattered-in process), and that the factor of two in front 
of the right hand side accounts for the sum over two photon 
polarizations (see Eq. |6]l. In the second equality, we have 
used the expression of the differential cross section given in 
Eq. (|53] |. where appears the squared amplitude |" of the in- 
verse Compton process. We further need to average to above 
expression, as previously done for the out part, over the solid 
angle dO. ~ d^d^, and the final expression for the scattered-in 
term is thus: 



fd^k in, . 



(17) 



As one can see, Eqs. ( fTOl i and ( fT6b are similar to Eqs. ( IB2l l and 
(|B3b . While the scattered-out processes are strictly equivalent, 
a slight difference appears in the scattered-in, coming from the 
squared amplitude (|rW |^ versus |5W |^) and the energy trans- 
fer (f versus t). Nevertheless, this difference is completely 
removed iii the relativistic Umit, since we have in this case 
rW p = |fW P and t = t. Hence the two methods are indeed 
equivalent. 



D. Redistribution function for the rSZ effect 



In Sect. 1111 Dl we will derive the full expression of the dif- 
ferential inverse Compton cross section for relativistic elec- 
trons and soft photons, and also provide some results in the 
corresponding limit a — Ej^/Ep (see Eq.l53Tl. Anticipat- 
ing these results, and using Eqs. ( 116117b . this implies that the 
angular average of I™{Ei^) can actually be written at leading 
order in a as: 

^ ^jdlj dpMEp,x) 

X j d^i' j d^i'F{^,ii,ii')f^{tEk), (18) 

where we define: 



(27t)3 16 (1-Pa')2 
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(19) 



K is defined in Eq. ( |49l ). The function f is meant 

to be related, while not strictly, to the product p{p) -qlii.ii') 
(where q{fj,ij') is in fact the squared matrix amplitude) which 
was displayed in Wright's paper However, it is different 
from Wright's expression, notably because this latter is de- 
rived in the electron rest frame, but they have been shown to 
be consistent ll2lll . 

Moreover, we see that in Eq. (fTSl l. we cannot factorize out 
the optical depth T as defined by Eq. ( fT2] i: it is mixed with 
the definition of T , and it cannot be taken out of the integral 
over the angles. Anyway, by re-expressing the photon zenith 
angles /j and fj' in the electron rest frame, it is possible to 



derive a much more suitable expression for numerical compu- 
tation ll2lll . consistent with the results obtained by Wright. 
Anyway, the trick one can use here to recover a formulation of 
the in process similar to the classical radiative transfer method 
used by Wright is the same as employed for the out process. 
Indeed, we can still describe the electron phase space density 
as fe{Ep,x) = ne{x)fe{Ep), where fe is normalized to unity, 
and utilize the expression of the non-relativistic optical depth 
Tnr given in Eq. ( fT6] l. Armed with this, we can express in 
contribution a la Wright as: 

= 2 Tnr Jdpjdijl dlj'9{^,IJ,^t')I^{tEk) (20) 



where we have defined: 
?(P,A',A'') 



Ot 



feiEp) . 



(21) 



These are some expressions complementary to those derived 
by in Wright'. The function f has to be interpreted as 
the photon frequency redistribution function for a given nor- 
malized energy distribution of relativistic electrons such that 

!d^p/{2%ffe{Ep) = 1. 

To summarize, the Boltzmann approach allows to derive con- 
tributions for the scattered out and /« processes which are sim- 
ilar to the radiative transfer method when regarding the rela- 
tivistic SZ effect. These contributions are recalled here: 



= 2Tn 



2^Xr.,^{Ek) 

jdpjdf,jdf,'9{^,ijMm) 



(22) 



We recall that the factors of two appearing above explicitly ac- 
count for the sum over the two polarization states of the CMB 
photons, so that ly must be defined for one degree of freedom 
only. Note that for numerical computations, the photon fre- 
quency redistribution function f is more suitably expressed 
when photon zenith angles are re-expressed in the electron 
rest frame, as shown by ll2lll . to which we refer the reader for 
computation purposes. 



E. Experimental signature of the deviation 

Finally, since an experiment has a finite angular resolution, 
the relevant quantity to use when doing predictions for spe- 
cific instruments is the average value of AIy{Eic) over the solid 
angle AHres carried by the experimental resolution, that is: 



d^res A/y(£'i:) 



(23) 



where the dependence on Hres of the intensity is hidden in the 
spatial dependence of the electron population characterized by 
fi,{p,x). More precisely, if the electron density is spherical in 
the cluster, then it only depends on the distance r = \x — Xcc | 
to the cluster center. This radius r can be related to the angle 
x|/res scanning the resolution range through the relation: 



r = yP + d^ -2c//cos\|/res , 

where d is the distance from the observer to the cluster center, 
and / is the distance as measured along the line of sight. 
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III. (INVERSE) COMPTON SCATTERING 

In this section, we rederive the famous resuhs for the squared 
matrix ampHtude and the phase space factors associated with 
the electron-photon elastic scattering process in the special 
kinematics of inverse Compton scattering. The number of 
references associated with the Compton scattering is obvi- 
ousl y very larg e, but one can find some more details in 
e.g.MlBlIi 

A. Kinematics 

The kinematics of the process we are interested in is the scat- 
tering of an electron of four-momentum p off a photon of four- 
momentum k, resulting in final states p' and k' for the electron 
and the photon, respectively. We deal with an elastic collision, 
in which the conservation of four-momentum holds, so that: 



p + k^p' + k' . 



(24) 



We can now derive almost everything we need from these ba- 
sic equations of classical collision theory and from the very 
simple geometry of the problem ll20ll . Note that the coming 
results will hold in any Lorentz frame, unless we specify it. 



1. The scattered photon energy 

The first quantity that we can calculate is the energy E/^i of the 
outgoing photon after the collision. Starting from Eq. ( l24b . 

— \/nv- + p'^, and p' — p + k — k' 



and substituting Epi 
find: 



m 



p 

l\2 



[p + k-k 

^Ey 

with t 



{Ep+Ek-Ek'f 
{Ep+Ek-Ek'f 



we 



(25) 



tE^ 



where P ; 



(l-p^') + a(l 
; p jEp is the electron velocity, where 



-A) 



a ; 



(26) 



which will be a quantity very useful to study the different as- 
trophysical regimes; and where we have defined the cosine 
of the angle 9 (9') between the incident electron and the initial 
(outgoing) photon by p (jj', respectively); and that of the angle 
between the initial and final photon directions by A: 



p = cos 9 = 
p = cos 9' 
A = cos© = 



p - k 

\p\-\k'\ 
k-k' 



(27) 



We can already have a look to the limit value of the scattered 
energy when the incoming electron energy is much larger than 
that of the incoming photon Ep Ek <X ^ Q, which is typ- 
ically the case in any collision between a CMB photon and an 
electron (indeed, Ei^^ <^ m): 



Ek ' (i-Pa'O 



(28) 



We see that, as p ^ 1 , /i^/ is maximized for frontal collisions 
scattering a photon of initial direction p~ — I to the opposite 
direction p' ~ 1. In this case, Ei^i — » A'f'Ei,, obtained after a 
Taylor expansion of p ~ 1 — 1 /(2'y^) up to the second order in 



2. Geometry 

It will be found convenient to express A as a function of the 
other cosines p and p' . Without loss of generality, we can 
choose the incident electron direction along the z axis of any 
frame Oxyz. Therefore, we can define the coordinates of the 
quadri-momenta as follows: 



p = £p(l,0,0,p) 

k ~ £'j;(l,sin9cos(|),sin(|)sin9,cos9) 
k' = £■,(./( 1, sin 9' cos (|)', sin (|)' sin 9', cos 9') , 



(29) 



where we have used p = \p\/Ep and the fact that = Ek for 
photons. The coordinates of p' are readily deduced from the 
above expressions by invoking energy-momentum conserva- 
tion. With these coordinates, we recover the previous defi- 
nitions of p and p' (see Eq. |27l ), and we can now express A, 
the cosine of the angle between the incoming and outgoing 
photons, as: 



A = cos© = 
where ^ = (|) — (j)' 



pp 



' + >A~^cos 



(30) 



Armed with these definitions and quantities, we can now sur- 
vey the full relativistic phase space as well as the squared am- 
pUtude of the process. 



B. Phase space 

The first relevant quantity when dealing with the probabil- 
ity of any quantum scattering process is the Lorentz invariant 
phase space c/LIPS, namely the possible final states to which 
the process can lead. With the same conventions as previously, 
we can write: 



c/LIPS 



d^p' d^k' {2Tlfh\p + k-p' -k') 



(27t)3 (27t)3 AEp,Ey 
By performing the integral over d^p ', we get: 



c/LIPS = 



d^k' b{Ep + Ek-Er~Ey) 



4(27t) 



Ef-Ek 



(31) 



(32) 
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where, accounting for the momentum conservation 
E„i > Ey, we define: 



Er = \Jm^ + {p + k-k'Y 



(33) 



El + El + El + lEp^iEkP - E^,^') - 2EkEk,A 



where the last line is obtained after having expressed /j as 
a function of t. Therefore, we see that the argument in the 
8 Dirac function of Eq.dSTTi is a non trivial function of Ej^i, 
f{Ei^i) = Ep +Ek — Er —Ei-i. By means of the properties of 
the Dirac functions, we can write: 

where El_ are the zeros of /, and /' stands for the derivative 
with respect to Ej^t. There is only one zero, Ey — E^ — t £4, 
given by Eq. dZSl l, so that we find: 



with E, 







(34) 



= Ep^J\+2a^{^i-t^i') + a?{\-2t^ + t^) 
= £p|l+a(l-f)|. 

Therefore, the Dirac function becomes: 



Epl-Sl 



(35) 



where « = a{t - A) - ^^i' + E^ /Ep 
= (l-p^') + a(l-A). 

We can go ahead with the derivation of the phase space factor, 
which can now be written as: 



dhlPS 



dD.ki k'^dk' E^^b{El - El,, ) 

2— -^r-TT^ — (36) 



E^, 4(27t)2 ErEk, Ep |« 

d^iii Ek t 



A{2%fEp |«| 



4(27t)2 



The second line is obtained after integration over k''^dk' 
ElidEki. 

As previously, it is interesting to consider the limit a (see 
Eq. |26] |. which is valid in the case of CMB photons scattered 
by hot or relativistic electrons {Ep ^ E^). In this limit, t 
(1 - P^)/(l - p^') and we finally get: 



E,,..p' 



dUFS ^ ^ (37) 



4(27t)2 (l-pA(' 



.^2 



This result describes the whole phase space function for elec- 
tron interactions with CMB photons. It can be expressed in 



terms of dt or dt, (where ^ = 1 /t), if we make the change of 
variable f <-> /j' or equivalently t ^ /j. With the definition of t 
given in Eq. dZSl l. and in the limit a ^ we have: 



, _ (1-Pa^) dt 

^ ^ P f2 



(38) 



dp 



dt 



Thus, in terms of f or ^ , the phase space is merely: 



J JT Tr,o a^o di!f' a , , 

du dLiPb > y-^ -^du dt , 

E,.p' 4(27t)2p 

a^O dd}' at . 

' 



(39) 



or / dpdUPS 

JEj., .p ' 



a^O 



4(27t)2 p 

diSf' a , d^ 
-dp 



4(27t)2p ^2 ' 

-,0 di^' a , di^ 
~^ 4(2^ p ^ ^ 



C. Squared matrix amplitude 

Only two Feynman diagrams have to be considered for the 
(inverse) Compton diffusion at leading order The amplitude 
of the process p + k ^ p' + k' is to be found in any textbook 
introducing quantum field theory (see e.g. lfl9ll '). and can be 
written as: 



-ie^Uf{p') Si Ui(p) , 



(40) 



with = ^-—J-^i^^i^ 

where / and / flag the initial and final states, p and k (prime) 
are the four-momenta of the incoming (outgoing) electron and 
photon, respectively, m is the electron mass and M,yj is the 
electron four-spinor, while e,yy is the space-like photon four- 
polarization. The Mandelstam variables are: 



s = [p + kf = n? + 2p.k = n? + 2EpEk{l-^p) 

--2p.k' ^ n? - 2EpEk'il - Pa(') 



(41) 



If we define = s — nP'= 2p.k, Du = u — rrP' = —2p.k', = 
p + k and qi, ^ p — k' , then the squared value of the amplitude 
in Eq. ( l40b . once averaged over the incoming electron spin, is 
given by: 



spin, pol 



2^2D,D„ 



Tr {d]Ti + Dll2 + D,D,'i^) 



with Ti = [P' + m)f,i{/iu+m)ff{P + m)if{li,,+m)fi 
+ (e,-^e/&& qs^qu) ■ 
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After some Dirac algebra calculation, and after averaging 
over the photon polarization states, one can find the very 
well known result for the squared Compton amplitude (see 
e.g. M): 



\M I 



+111* 



P-k ^ p-k' 
p -k' p-k 



(43) 



1 



1 



p-k p-k' 



1 



1 



p-k p-k' 



Using the conventions displayed in Sect. IIII A 21 we have the 
four-products: 



p-k = myEk 



and 



p-k' = mjE,, (I-Pa/') 
so Eq. ( |43T l finally reduces to: 

(1-A) 



(44) 



(45) 



= 2e*{l 



1 - 



y2(l-pA')(l-pA'0 

£2(1 -A)2 



y2(i-Pa'0[(i-Pa'') + «(i-a)] 

2 



} 



2/ <!2 



y2 \ y2 



1+ag 



(46) 



where we define: 



g ^ (l-A)/(l-p^') 
and f = y^(l -P^) . 



(47) 



Y = Ep/m is the usual Lorentz factor for the incoming elec- 
tron, and e = E/^/m = ay is the incoming photon energy in 
units of the electron mass. The second equality given in 
Eq. ( |46] l is very convenient because the result of the limit 
Ep ^ E/t <^ a ^ appears explicitly. Within this limit, the 
squared amplitude is merely: 



1 2 a^O 



2e^ 



y2 



(48) 



Should we implement the non-relativistic limit P ^ 0, we 
would immediately recover the well known angular depen- 
dence of the Thomson scattering, that is |rW |2 oc 1 + a2. Any- 
way, the above expression depends on the azimuthal angle 
through A (see Eq. [30] l, and it is interesting to perform the az- 
imuthal integral for deeper insights. It is clear that only the 
even powers of cos^ will contribute to this integral. We find: 



with K 



2/ I^2-2K{1 -iJiLt') 



(49) 



[(1-A<A'')'+2(1-A'')(1-A'") 



-1 



From this expression, it is easy to calculate the limit P — * 0, 
that is the expression corresponding to the scattering of a pho- 
ton of energy Ek off an electron at rest. Indeed, within 
this limit, K ^ \, and the squared amplitude becomes: 



(|^|2)^ ^ 2e^ 



This is the famous non-relativistic limit derived by Chan- 
drasekhar in his textbook ||20|] and used in most of treatments 
of SZ computation, and in particular in (by Lorentz boost- 
ing the photons up and down with respect to the electron rest 
frame). This expression can not be used such as it is writ- 
ten above in relativistic SZ calculations, but it turns out that 
a form very similar can be derived when expressing the pho- 
ton zenith angles /j and /j' in the electron rest frame, that is 
A- = (P -W)4.1 - Pa-o) andA<' = (P -A'J,)/!! - pA-o), as nicely 
detailed in 112111. 



D. Differential (inverse) Compton cross section 

The (inverse) Compton differential cross section is defined as 
usual for two body elastic scattering {p,k —> p' ,k'): 



(2^ 
Ap-k 



n 

^X=p'.,k' 



{2%f2Ex 



■b\p + k-p' -k')\M\^($l) 



where \ already includes the sum and average over polar- 
izations, and where the 4-product p-k = EpEk%e\ in terms of 
the energies of the incoming particles and their associated rel- 
ative velocity. After integration over p ' and Ej^t, by virtue of 
the energy-momentum conservation, the angular differential 
cross section reads: 



d<5 



1 



c/LIPS ISW 



dQ.1,1 4EkEp^rel Jp',Ei^, 

After substituting prei = (1 - p^*) (see Eq.l44b: 



da 



1 



647t2F2(l_p^)2 



(52) 



(53) 



3ot 



a^O 



167IY^(1-Pa')^ 2e4 
3ot 

167172 

2 



+K 



^{2-2^(1-..') + 
[(l-^A'r + ^(l-A'')(l-A'")]} 



where we have used the standard definition of the Thomson 
cross section Oj = Snr^/S, with the classical electron radius 
ro = e^/{4-Km) = ttfj/m. The first two lines give the exact 
and general expression of the (inverse) Compton scattering 
differential cross section. The full relativistic expressions of 
t and \!M p have been derived in Eqs. ( [25] ) and ( |46l ), respec- 
tively. The third line shows the relativistic limit at leading 
order, when a ^ E^/Ep ^ 0, including the average over the 
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photon azimuthal angle, according to Eq. ( |49] l where K is also 
defined. If we further integrate' the full differential cross sec- 
tion over dD-k', we find (see e.g. ll27[l28ll '): 

o(£p,£.) = ^|fc^^ln(l+2x) (54) 

(l-f2x)^+(l+X)'(l+2x)-x' l 
X2(1+2X)2 / 

> Gt , 

where % = ay^(l — P^). We have featured the limit when 
ay^ = yEk/m 0. It makes sense to recover the Thomson 
cross section, since if we transpose the process in the electron 
rest frame, the photon would have an energy of yE^, still much 
lower than m as long as y < 10**, which corresponds to the 
classical Thomson regime. While not describing a Thomson 
scattering at all, the corresponding integrated cross section is 
recovered as long as yE^. <C m: this is a mere consequence 
of special relativity, and very well know for decades in astro- 
physics imii. 

Another useful quantity to derive is the integral over the in- 
coming and outgoing photons angles of the Taylor expansion 
of pieit/o/c/n^./ in the limit a ^ 0. Up to the second order in 
a, we get: 

(55) 

where we have expressed the cosine A as a function of fi, /j' 
and (|)' according to Eq. dZTl ) before performing the integral. 



IV. CONCLUSION 

We have shown that the approach that Wright developed in 
Ref. ^ was equivalent to using a Boltzmann-like equation 
in the single scattering approximation, already very well de- 
scribed in the literature (e.g. ifiol [Til [l3i [Till ). The same result 
has been recovered in a recent subsequent work by Nozawa 
and Kohyama 1I21I1 . Here, we have established the relativis- 
tic expressions valid in the CMB frame to employ within 
Wright's formalism, summarized in Eqs. ( l22l i. They differ 
from Wright's expressions which are instead derived in the 
electron rest frame 1I21I1 , and which are probably more suited 
for numerical computations. Finally, we have recovered that 



We temporarily slightly change our convention which defines the cosine A 
with respect to /j and /j' — see Eq. i27\ . We now set the main axis of the 
problem along k ' instead of p. so that it is now fA which depends on A and iJ 
and so that we can perform the integral over d/J without accounting for A. 
Such a trick is aimed at comparing our results with the standard calculation 
of the integrated Compton cross section. 



the SZ shift at energy Ek due to a relativistic population of 
electrons is given by an integral over a function which is still 
proportional to Iy{Eii) — Iy{t Ek)lt^, where /y" is the intensity 
of the pure black-body CMB spectrum and t = E^./ /Ej^ is the 
ratio of scattered-to-focused frequencies. 
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APPENDIX A: RECALL OF THE RADIATIVE TRANSFER 
APPROACH 

In this appendix, we would like to demonstrate that the ra- 
diative transfer method used by Wright in ^ can be deduced 
from the Boltzmann approach presented in this paper, but by 
starting from the radiative approach instead. This will provide 
more insights on the correct frequency redistribution function 
to be used for the CMB photons. 

We are looking for a redistribution function that expresses 
the probability P{t) for a CMB photon to have its energy Ek 
shifted to E^j by a factor of f = E^i /E^ — we could equiva- 
lently deal with P{s) with s = ln(f ) — in order to be able to 
write (see e.g. U^) the modified intensity as: 

Iy{Ek) = Jdt P{t) X I^itEk) . (Al) 

The integral is performed from the minimal value of the scat- 
tered energy fmi,, = 1 / (4ymax — 1 ) up to its maximal value 
tmax ~ (4l^ax ^ 1)' allowed by the kinematics. 
According to the radiative transfer method, such a probabil- 
ity depends on the electron optical depth T and on the fre- 
quency redistribution function Pi (?) for a single scattering. In 
the limit of single scattering, it is merely given by: 

P{t) = e-'{5{t-l)+zPi{t)} 

^ (l-T)8(f-l)+TPi(f) + 0(T2), (A2) 

where, in the second line, we have taken the limit of very small 
optical depth, which is the case in the problem of interest here, 
but also more generally for electrons in clusters. 
If we report the previous equation in Eq. ( lAlk we find: 

AlyiEk) =^Jdt (Pi (f) - 5(f - 1)) X I^itEk) . (A3) 

Two terms appear in the sum in the integral, which are readily 
interpreted as the positive contribution of CMB photons scat- 
tered in from energy E^i = tE^ to energy E^ for the former. 
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and as the negative contribution of photons which have been 
shifted out from energy Ej^^ to other energies. More explicitly, 
those contributions are given by: 

/7'(£i)=^x/0(£,), 

and: 

=^jdtP, (f ) X /0(f£,) . (A5) 

If we scrutinize Eqs. (|9]l and ( |20l i. we can find some convinc- 
ing similarities, while not straightforwardly. With the defi- 
nition of T given in Eq. ( fT2] i. the former is perfectly consis- 
tent with Eq. ( IA4b . Nevertheless, it is a bit more tricky for 
the latter Indeed, given our definition, because T already in- 
cludes the integral over f, it would have to be considered as 
an operator acting on /y. We have actually determined the full 
relativistic formula, summarized in Eqs. (l22l i. 



APPENDIX B: SCATTERED IN AND OUT PHOTONS 
TREATED SEPARATELY 

1. Definitions 

Let us start with some basic definitions. The energy density 
of the black-body spectrum, of average temperature 7b, before 
traveling through the cluster is given by: 

p« = | d'kE,f^\Ek), 

with d^k — E^ dEii dQ.k and 

f^^^''^ " eEkiT, _ 1 ■ 

From this equations, we deduce that the associated intensity 
is: 

/0(£,) = Elf^iEk) (Bl) 

gEklTo - 1 ■ 



2. Intensity associated witli tiie relativistic effect 

When interacting with the relativistic electrons, some CMB 
photons initially at energy E^ can be scattered out at an en- 
ergy Ej^i. This effect generates a deficit in the number of CMB 
photons that is expected at an energy Ei^. However, other CMB 
photons with an energy Ej^i ^ Ek will be scattered in the clus- 
ter to a final energy Ej^. The photons associated with the for- 
mer process (responsible for the deficit) will be referred to 
as "scattered out" photons (and will be denoted hereafter "s- 
out") while those associated with the second process will be 
referred to as "scattered in photons" (and will be denoted "s- 
in"). Let us now examine these two processes separately. 



a. Intensity associated with the photons scattered out 

After traveling through the cluster and interacting with the 
population of relativistic electrons, a CMB photon (with an 
initial energy E^) can be shifted to a new energy Ey ^ Ek, ac- 
cording to the reaction k+p^kf +p'. The CMB energy spec- 
trum is modified accordingly. For an observer who aims at 
measuring the resulting intensity of the CMB spectrum, such 
interactions with a population of relativistic electrons trans- 
lates into a deficit in the number of photons at an energy E^ 
with respect to the black-body distribution. 
If we disregard the thermal and kinetic SZ effect for a mo- 
ment, the number of photons which remain at this energy E^ 
after traveling through the cluster is the number of photons 
which had this energy before entering the cluster minus the 
number of photons «y^°"' which have changed of energy after 
traveling through the cluster. 

This number density «y^°"' of scattered photons is basically 

the convolution over the 3 -momentum d^k between the ini- 
tial black-body occupation number fy and the interaction rate 
with the electrons (that is Fy = «j.apiei if Opiei is independent 
of the ingoing electron energy Ep). For relativistic electrons, 
the dependence of the cross section on the energy Ep is not 
trivial and the interaction rate is rather given by: 

f d^P 

Ty{Ek) = J fe{Ep) [oPrel]pi:^p'A:' 

where 

with Prei the relative velocity and where the Lorentz invariant 
phase space c/LlPS = dLWSpi y is defined in Eq. dSTT i. and 
further developed in Eq. ( [36] l. |5Vf p = |5W 1^^^^/^,/ is given in 
Eq. ( |46] l. To derive the full distribution function of the photons 
which have been scattered out from energy Ek to energy Ek', 
we have to multiply (or weight) the black-body distribution by 
the interaction rate integrated over the time: 

/f = fy{Ek) j dtYyiEk) 

= f^{Ek) jdt j ^f,{Ep) [a^rel]kp^^p', 

which provides the fuU number density of photons scattered- 
out from any energy Ek to any Ek', that is the number density 
of photons having interacted: 

Now, we are interested in the specific density removed from a 
given energy Ek and shifted to any energy Ek', so we have to 
come back to the distribution function. Hence, the intensity 
associated with the scattered-out photons reads: 

I^-°-\Ek)^Elf;-°-\Ek), 
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where /^^""'(^i) represents the distribution function associ- 
ated with the photons which have experienced elastic scat- 
terings with electrons in the cluster Since an observer will 
detect a modification of the photon distribution along the line- 
of-sight, the deficit in the CMB intensity at an energy Ek that 
could be observed is merely obtained by using the geodesic 
relation dt ^ dl / c {— dl in natural units): 



shifting photons from energy Ej^t to energy Ek reads: 

dTy{Ei,,,Ek) f d^p dOpk'^p'k 

^ riJe{Ep)Prel 



d^k 



(2%) 



d^k 



dl 



d^p 

d^p 



fe{Ep] 



rfLlPS^^/0(£, 



To derive the distribution function associated with all photons 
scattered from energy E^i to energy E^, we need to convolve 
the above differential rate with the initial black-body function 
and integrate over time and over all energies E^i: 



647^2 (I-Pa') Ej, 



where the Lorentz invariant phase space has been made ex- 
plicit. We have finally: 



\Ek) 



dl 



d^p 
j2Kp 



fe{Ep) 



d^L 



327:2 (1 -^ 



£■2 



(B2) 

where we used the fact that /^(-Ei) = 2E^ /^{E^), making ex- 
plicit the factor of two coming from the sum over the photon 
polarization states, and where the squared amplitude |5Vf p = 
We exactly recover the scattered-out term as de- 



p,k^p' .k' 



rived in the Boltzmann-like formalism (see Eqs.[T0land[T2b. 
Of course, one will have to average over the solid angle of the 
incoming photons dflk ~ di^dfj. 



b. Intensity associated with the photons scattered in 

Let us now estimate the number of photons with an arbitrary 
energy Ek' which have actually been scattered to an energy 
Ek by a population of electrons of energy Ep. A first approach 
consists in doing as in the previous subsection, by defining the 
interaction rate of scattering shifting photons of energy Ek' to 
any energy Ek ^ Ek'; 



TyiEk') = 



d^p 



fe{Ep) [0^re\]pk'^p'k- 



^0/z7 N f jrlj, NO dOpk'^p'k 



(271 



3^7 1 



(27^)3 



d^k 



We have: 



dOpk'^p'k (271)4 . ^3^/ 



'AEkEp, 



where |5Vf p = \'M l^pk'^p'k- '^hat the phase space is 

more subtle than in the scattered-out term. We can actually 
perform the integral over d^k' and d^p', by means of the con- 
servation of energy-momentum. This will differ a bit from the 
scattered-out case, but we will apply the same treatment as in 
Sect. 1111 B1 to which we refer the reader for more details. We 
define: 



dLlPS 



This is not the common Lorentz invariant phase space for the 
process pk' p'k, which should instead feature integrals over 
p' and k, but this is a useful quantity in the present calculation. 
If we perform the integral over d^p\ we get: 



dlA?S 



d^' 



\6%^EkEr 



-r^Ep + Ek' - Er - Ek) 



with: 



Er = EpJ 1 + a''^ + a? + 2^{a'/j' ~a/j)-2aa'A , 



The number density of scattered-in photons from any energy 
Ek' to any energy Ek is therefore given, in terms of the associ- 
ated distribution function, by: 



d^k' 



^ fyiEk') 



f^{Ek') = f^{Ek') I dtTyiEk 



Nevertheless, this is the number density of photons scattered 
to any energy Ek ^ Ek', that is the total number density of 
photons having interacted, while we are only interested in 
the scattered photons from any energy Ek' to a given energy 
Ek- Therefore, we have to modify a bit the above picture, and 
use differential quantities. The differential interaction rate of 



where a = Ek/Ep, a! = Ek' /Ep, and the other angular vari- 
ables are already defined in Eq. ( |27] |. We can further transpose 
the Dirac 5 function according to: 



with: 



d{Ep+Ek' -Er-Ek) 



E^m'-E'k) 



Ep\'B\ 



El = lEk 



(l-pA'') + a(A-l) 
« = (1-Pa;') + oc(A-1) 
£0 ^ E,-{El)=Ep\l-a{l-t)\ 
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Hence, performing the integral over cf'k' , we can rewrite the 
phase space as: 



such that the scattered-in distribution function can be written 
as; 



647:2 El (1-p^) 



We can now derive the intensity of the scattered-in photons at 
energy fi^t; 

If we transform the time integral into a line of sight integral, 
we finally get: 



\Ek) 



dl 



d^p fdSly |5W|2 I^yim) 



327^2 £2 ,~(i„p^) 



(B3) 

where we have used Iy{tEk) = 2{tEkf' f!^{tEk) — making ex- 
plicit the factor of two coming from the sum over the pho- 
ton polarization states — and where the squared amplitude 



|rW |2 — |rW Of course, we will further need to per- 

form the averaging over the solid angle d£lk = di^d^i associ- 
ated with the detected photon. Nevertheless, it is already in- 
teresting to note that we recover an expression which is simi- 
lar to Eq. ( fTSI l, which is the result for the in process obtained 
in the Boltzmann-like formalism, but with slight differences. 
Indeed, in the Boltzmann approach, the squared amplitude is 



\'M \pi^^pH.t, and the expression of f = jE^ slightly differs 
from that of t. As regards the energy ratios, the difference 
comes from the denominators: there is a factor of a(l — A) 
in this of the former, and of — a(l — A) in that the latter. The 
physical interpretation of this slight difference is rather sim- 
ple. Indeed, a = Ej^jEp in both cases, but though Ef^ is the 
energy of the incoming photon in the former case, it is that 
of the outgoing photon in the latter case. Nevertheless, since 
a — ^ for CMB photons as scattered by relativistic electrons, 
this has no effect. Concerning the squared amplitudes, we 

can also verify that |rW |2 and \'M |2 are equivalent within this 
limit. There is therefore no difference between the scattered- 
in contribution as computed in the approach developed in this 
appendix and that computed within the Boltzmann-like for- 
malism in the limit a ^ 0, which fully applies for relativistic 
electrons. 
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